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Recently it has been proposed to construct quantum error-correcting codes that embed a finite- 
dimensional Hilbert space in the infinite-dimensional Hilbert space of a system described by contin- 
uous quantum variables [D. Gottesman et al., Phys. Rev. A 64, 012310 (2001)]. The main difficulty 
of this continuous variable encoding relies on the physical generation of the quantum codewords. 
We show that ponderomotive interaction sufhces to this end. As a matter of fact, this kind of 
interaction between a system and a meter causes a frequency change on the meter proportional to 
the position quadrature of the system. Then, a phase measurement of the meter leaves the system 
in an eigenstate of the stabilizer generators, provided that system and meter's initial states are 
suitably prepared. Here we show how to implement this interaction using trapped ions, and how the 
encoding can be performed on their motional degrees of freedom. The robustness of the codewords 
with respect to the various experimental imperfections is then analyzed. 

PACS numbers: 03.67.Pp, 32.80.Lg, 03.65.Ta 



Quantum information, as well as classical information, can be carried either by discrete or continuous variable (CV) 
systems. The latter have attracted an increasing interest during last years Q. In particular quantum error correction 
(QEC) techniques have been extended to this framework in order to allow CV quantum computation |2j,|3j]. However, 
CV QEC presents some difïerent aspects with respect to that concerning qubits. In fact, while it is reasonable to 
encode a qubit in a block of qubits using a QEC code protecting against a large error on a single qubit of the block, 
this is no more true in CV systems. In such a case, in fact, the most probable effect caused by decoherence is a 
small diffusion in the position and momentum of all the partides. This is just the kind of noise treated in where 
shift-resistant quantum codes for qudits are suitably extended to CV systems. The main drawback of such codes is the 
preparation of the encoded states which ideally are non-normalizable states and therefore can only be approximated 
introducing an intrinsic error probability. Travaglione and Milburn Q| have shown that it is possible to generate such 
approximate states by performing a sequence of operations similar to a quantum random walk algorithm 5]. More 
recently, we have proposed an all-optical scheme based on the cross-Kerr interaction to realize such a code [6j • In this 
paper we follow the suggestion of Ref. and we study a ponderomotive interaction 7] to embed a qubit in a CV 
quantum system. 

The paper is organized as follows. In Sec. U^we rapidly review some elements from Ref. and we turn from an 
ideal situation to a more realistic one. In Sec. lIIIl the physical implcmcntation of the CV encoding scheme is proposed. 
Sec. II VI is for conclusions. 



A single qubit living in a Hilbert space Ti. with basis {|0) , |1)} can be encoded into a single oscillator in such a 
way that the two resulting codewords |0), |1) provide protection against small diffusion errors in both position x and 
momentum p (the quantum operators obey the commutation rulc [x,p] = i so that x,p are dimensionless quantities). 
The two quantum codewords |0), |1) are the simultaneous eigenstates, with eigenvalue +1, of the displacement oper- 
ators D x {29) = er 2l9 P, D p (2ir0- 1 ) = e 2 ™ 9 ' 1 *, with 9 € R, which are also the stabilizer generators of the code 0. 
These codewords are therefore invariant under the shifts x — > x — 29 and p ^ p — 2ir9~ 1 . Up to a normalization factor 
they are given by 
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II. ENCODED STATES 
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i.e. they are a coherent superposition of infinitely squeezed states (position eigenstates and momentum eigenstates). 
Each of them is a comb-state both in x and in p with equally spaced spikes (29 in x and tt9~ 1 in p). The codewords 
|0),|1) are also eigenstates of the encoded bit-flip operator Z ~ D p (nQ~ l ). Equivalently one can also choose the 
codewords |±) = [|0) ± \l)]/y/2 which are the eigenstates of the encoded phase-fhp operator X = D x (9) and are given 
by: 
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Also these states are comb-like states both in x and in p, with equally spaced spikes (9 in x and 2ir9 1 in p) . The 
four codewords states are schematically displayed in Fig. ^ 
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FIG. 1: Ideal encoded states |0),|1) (Z eigenstates) and |+),|— ) [X eigenstates). On the left the structure of the spatial 
wavefunctions is displayed while on the right the structure of the momentum wavefunction is displayed. Each spike is ideally 
a Dirac-delta function. 



The recovery process is realized by measuring the stabilizer generators D x (29), D p (2ir9 1 ). The measurement of 
the X-generator D x (29) = p(modir9~ 1 ) reveals momentum shifts Ap which are correctable if \Ap\ < ir9~ 1 /2; in such 
a case the correction is made by shifting p so to become equal to the nearest múltiple of ir9~ 1 . In the same way, the 
measurement of the Z-generator D p (2tt9~ 1 ) = x(mod9) reveals position shifts which are correctable if \Ax\ < 9/2; in 
such a case the correction is made by shifting x so to coincide with the nearest múltiple of 9. 

Ref. |3| proposed the following recipe for the generation of the codeword states. 

1. Preparation of a particle in the p = eigenstate (i.e. completely delocalized in position). 

2. Coupling the particle to a meter (i.e. an oscillator, with ladder operators c, &) via the non linear interaction 
Hnl = gc^cx. This interaction modifies the frequency of the meter by Alo = gx so that, at time t — ir9~ 1 g~ 1 , 
the phase of the meter is shifted by A<f) = tt9~ 1 x. 

3. Reading out the phase of the meter A</> at a time í, i.e. measuring x(mod29). This measurement projects the 
initial state into a superposition of equally spaced delta function S(x — 29s + e) with s = 0, ±1, ... and eel. 

4. Applying a suitable transformation to obtain any desired encoded qubit state a|0) + ò|l). 

Ideally the codewords are non-normalizable states infinitely squeezed both in x and p, but in practice one can 
only generate states with finite squeezing, i.e. only approximate codewords: |0) ~ |0), |1) = D x (9)\0) ~ |1), |±) = 
[|0)±|l)]/A/± ~ |±) ( Aí± are normalization constants). For this reason, in order to estimate the quality of the encoding 
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scheme, together with the error probability in the recovery process due to the occurrence of an uncorrectable error, 
we have also to consider the intrinsic error probability due to the imperfections of the approximate codewords which 
can lead to an error even in the presence of a correctable error. Here, we propose a physical implementation of the 
ideal coding protocol of Ref. 3], based on the interaction between an atom and a radiation mode. In general, this 
scheme can be derived from the ideal one by replacing the initial p — state with a hnitely squeezed state, Hnl with 
a ponderomotive interaction, and the phase measurement with a homodyne measurement. 

III. PONDEROMOTIVE ENCODING 

Our proposal for a physical implementation of the CV encoding scheme is schematically depicted in Fig. 

harmònic trap 

i i 
i i 



', ; homodyne 
,' detector 




> X 

FIG. 2: Ponderomotive encoding. An ion is trapped by a harmònic potential within a high finesse cavity, where it interacts 
with a single mode. After a suitable interaction time, a homodyne measurement of the intracavity quadrature X — (c + c*)/a/2 
projects the state of the ion onto an approximate comb-like state, i.e. the approximate codeword |0). A conditional displacement 
can then be used to generate any correctable state a\0) + b\l) (see text). 

We consider a single two-level ion (of mass M and infernal transition frequency lüq — 2ttc/Xq) trapped inside a 
high finesse cavity by a harmònic potential, with trapping frequency ui a along the direction of the cavity axis x (see 
Ref. for experimental schemes of this kind). We assume that the trapping potential in the other two directions 
is steep enough to freeze the motion along y and z. The ion interacts with a single cavity mode with annihilation 
operator c and frequency oj c — 2irc/\ c — ck c . If the harmònic potential minimum is set halfway in between a node 
and an antinode of the stationary field of the cavity mode, the ion-cavity Hamiltonian can be written as (assuming, 
as usual in the optical domain, the rotating wave and dipole approximation) 

H = hiü a a i à + tko c c'c + hujQ<j z + hgo (à'c + &à\ cos(fc c x + 7r/4). (5) 

In this Hamiltonian à z — (à^à — <7<7^)/2, <r, are the atòmic spin-1/2 operators associated with the two infernal 
levels whose transition is quasi-resonant with the optical cavity mode, x = ^h/(2Mcu a )(à + aJ) is the atòmic center- 
of-mass (CM) position operator, c, are the cavity mode annihilation and creation operators, and go is the atom-ficld 
coupling constant. 

A ponderomotive interaction between the atom and the cavity mode is obtained in the dispersive limit in which 
the cavity mode is highly (red) detuned from the atòmic transition. In this limit, the upper atòmic level can be 
adiabatically eliminated and also the spontaneous emission from it can be neglected. In such a condition the atom 
always remains in its ground state and the resulting ponderomotive Hamiltonian (in interaction picture with respect 
to Ho = HlücC^c), is 

H = huj a à^à + H^-c^c cos 2 (k c x + n/4), (6) 
ò 

where S — ujq — lo c . We then make a second assumption, the Lamb-Dicke limit, which amounts to assume that the 
ion CM is well localized with respect to a wavelength of the cavity mode A c . This assumption allows us to linearize 
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the resulting optical potential in the far-off resonance regime. We thus end up with the simple Hamiltonian [ïo| 



H = huj a à^ à — h^zc) c — hgc^ c (à + à^) , (7) 
2o 



where g = (<?o/£)£ and £ = k c ^/h/2Muj a is the Lamb-Dicke parameter. 
The time evolution operator takes the forní 0, 0] 

Ü(t) — e ik2 ( £Í £ Y ( T - sinT ) e k£Í È (n à ^ -V* à) ç-irà 1 à ^ 

where r = u> a t, k = g/uj a , r\ = 1 — e~ ÍT . In writing down Eq.JH|) we have disregarded the free radiation field evolution 
term. 

We now consider as initial condition of the ion-cavity system the tensor product of a coherent state of amplitude 
a for the cavity mode and a position squeezed state of amplitude /3 and squeezing parameter e = r exp {2icf>} for the 
ion center of mass (see Ref. for the preparation of motional squeezed states of the ion CM) 

I*(0)) = |i9,e) a ®|a) c , \p,e) a = D(J3)S(e)\Q) a , (9) 

where D (/3) is the displacement operator and S (e) is the squeezing operator |14| . In order to compute the time- 
evolved state (r)) = U (r) \^ (0)), we expand the optical mode coherent state \a) c in number states \n) c and use 
the Baker-Campbell-Hausdorff expansion [Ï3. After some àlgebra, we obtain 



|tf (t)) = ^-^- exp 
Vcosnr 



r/3* 



oo 

X 



^= exp[r/3*e- lr (tf - nfcr;*) - (r/2)e- 2ÍT (a f - n/by*) 2 + ik 2 n 2 (r - sinr)] 

x |/3e" lr + nkrj) a <g> |n) c , (10) 

where |/3e _lr + nfcr?) are coherent states of the ion CM motion and T = exp {2ï</>} tanhr. Inserting the identity 
decomposition in the basis of ion CM coherent states, I = J d 2 j \j) a (j\ /tt in Eq. l(TÜ|) . we obtain 

|* (r)) =Aj2^n — e c - | 7 ) <8> |n) c , (11) 
rrí J * 



n=0 



where we have defined 



with 



A= (coshr)~ 1/2 e~ 2"" , (12) 

C„, 7 = [r(/3* + C) + /3 + Cn]e- lT 7* - (I72) e - 2 *V 2 " |7| 2 A (14) 

C n = nfc(e ÏT - 1). (15) 

At the (scalcd) time r we measure the intracavity quadrature X — (c + c^)/\/2 obtaining the result X 0. As 
a consequence, the cavity mode is projected onto the corresponding quadrature eigenstate \X), while the ion CM 
motion is disentangled from the cavity mode and it is projected onto the state with wave-function 

N °° 

*(*) = ^w==^ ^E ( X \n) c B n e^ , (16) 

v n—O 

where Nx, r is a normalization constant, and x now stands for the dimensionless ion CM position, i.e. x — > x^J Mu) a /h. 
Furthermore 

P»W = ~(r - e^J-^ír + e 2i > 2 + [T(/T + C) + P + C»][r(/3* + C) +P + C» - 2V2V T x]} , (17) 
with (X|n) c = tt- 1 / 4 (2™n!r 1/2 H n (X)e- x '^ 2 , H n (X) being the n*' 1 Hermite polynomial. 
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A. Homodyning with a zero outcome 



We have now to determine the conditions under which the general ion CM conditional state of Eq. I|16fl becomes 
a comb-like state which can be taken as approximate codeword state. It is possible to see that one needs to choose 
k = 1/2, that the homodyne measurement has to be performed at the appropriate time t = tt, and that initially the 
ion has to be displaced and squeezed in position, which is achieved if we take (3, e real and positive. To see this we 
consider a particular example of homodyne measurement outcome, which makes things easier to see, i.e. X = 0. In 
such a case in fact, the property i?2m(0) = (— l) m 2" l (2m — 1)!!, i/ 2m +i(0) = (m = 0, 1, ...) allows to get from lfTÏÏ|) 
the following expression for the ion CM wave-function 



tp(x) = $(aO| fe=1/2 . T=7r;/3)e >o ; X=0 = N Vrntt m (x) , (18) 

where 



m=0 



-a 2 /2_^ 



(2m)! 



(19) 



e 



r/2 



n m (x) = exp{-i e 2r [x - V2(2m - f3)] 2 } , (20) 
V 7 !" 

(we have also taken a real and positive) and N = [P(X = 0)] with 

oo 

V{X = Q) = — V VmVm ,e- 2(? ^ m - m '^ (21) 
Wtt í — ' 

m,m — 

being the probability density corresponding to the outcome X = of the homodyne measurement. The state of 
Eq. i(T%|) represents, in fact, a superposition of position squeezed states centered in x — \f2(2m — (3) , (m = 0, 1, 2...), 
which is however also a superposition of squeezed states in the momentum space p, as one can verify by computing 
the Fourier transform 

V>(p) = / <p(x)e~ lpx dx 



2ir 

exp[-(r + e- 2 V)/2] V v m e-^ 2m - « (22) 



(see Figs. Eíc), Efe), Eïg))- Therefore, also in the present scheme, the homodyne measurement of the cavity mode 
conditionally generates the desired comb-like state, which we can take as the approximate codeword state |0). This 
state depends upon the three dimensionless parameters a, (3, and r. The parameter r is obviously responsible for 
squeezing in x and increases the number of spikes in p (as we can see from Eqs. I|20|) . and (I22() '). while a causes 
squeezing in p and increases the number of spikes in x. This fact can be seen if we consider the probability to have 
the m th spike ft m (x) in the wave-function <p(x), P{m) — ^m/X)m=o Vm - ^ ac ^ one can see that the mean value fh 
and its Standard deviation Am are approximately given by 

m~a 2 /2, Am~a/V2, (23) 

implying that the number of spikes in x grows linearly_with a. These arguments show that the generated approximate 
comb-like state approaches the ideal codeword state |0) of Sec. [H]for a, r — > oo and that we have to take these latter 
parameters as large as possible. By varying the parameter j3 instead, one simply shifts the state along x. 

Therefore, by comparing with the ideal codeword state of Eq. we have that the state of Eq. i|18|) approximates 

it (except for an unimportant shift by (3) with 9 — y2. The approximate codeword |1) = D x (— v2)|0) is generated 
by a position shift equivalent to a change of the parameter /3, A/3 = 1, and which can be again realized by applying a 
suitable elèctric field (or làser) pulse (see Ref. 13] ). After seeing how the two basis codeword states are generated, let 

us now see how to generate a genèric superposition of the two, a|0) These superpositions can be generated using 

conditional displaccmcnt schemcs analogous to those used, for example, in the manipulation of quantum states of 
trapped ions |l3| and which exploit the coupling of a motional degree of freedom with an infernal transition of the ion. 

Schematically these schemes proceed as follows. The atom is prepared in the tensor product state |0) ® [a|<?) + b\e)], 
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FIG. 3: Plot of the wave-functions in space and momentum coordinates of the approximate codeword states |0), |1), |+) and 
|— ) for a — 1.8, r = 1.5 and (3 — 0. (a): spatial wave-functions ipo of |0) (sòlid line) and ipi of |1) (dashed line) vs x/\^2; (b): 
spatial wave-functions ip+ of |+) (sòlid line) and <p- of |— } (dashed line) vs x/v2; (c): momentum probability distributions 
\ïpo\ 2 of |0) (sòlid line) and l^i] 2 of |1) (dashed line) vs V2p; (d): momentum probability distributions \Í>+\ 2 of +) (sòlid line) 
and \ip- 2 of |— } (dashed line) vs \/2p; (e): imaginary part of the momentum wave-functions Imtpo of jO) (sòlid line) and Irm/>i 
of |1) (dashed line) vs V^p; (f): imaginary part of the momentum wave-functions lmtp + of |+) (sòlid line) and Imt/)_ of |— ) 
(dashed line) vs V2p; (g): real part of the momentum wave-functions Re^o of |0) (sòlid line) and Re^i of |1) (dashed line) vs 
\[2p\ (h): real part of the momentum wave-functions Reip + of |+) (sòlid line) and Retp- of |— } (dashed line) vs \[2p. 
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where |e) and \g) are two ground state sublevels. Then a làser pulse which is only coupled to |e) is applied to the atom 
and its intensity is tuned so to give exactly a position shift x — > x — \Í2. In this way the state of the atom becomes 
a\g) <g> |0) + b\e) ® |1). Then a rf pulse resonant with the e — > g transition and transforming \e) — > (|e) + \g))/y/2 and 
\g) -* l e ))/\/2 is applied, so that the state of the atom becomes [\g) <g> (o|0) + 6|1)) + |e> ® (6|1) - a|0))]/\/2. When 

the internal state of the atom is measured and it is found equal to \g), the atòmic motional state is conditionally 
generated in the desired encoded superposition a|0) + 6|1). Particular examples of these superpositions are the 
eigenstates of the encoded phase-flip operator X, whose approximate versions are given by |±) = [|0) ± |l)]/A/±, 
where J\f± — 2[1 ± (0 11)] (|0) and |1) are not orthogonal in general also in this scheme). It is interesting to notice that 
the corresponding wave-functions of these latter codeword states |±) in the momentum space are given by 

oo 

1> ± (p) = {p\±) = n-^NN^ e 'V2/3 P -(-+e- 2 V)/2 (1 ± e V2 P) u m e~ i2 ^ mp (24) 

showing that \tp+ (p)\ 2 has spikes at p n = 7r(2n)/\/2 while \ip_ (p)\ 2 has spikes at p n = n(2n + \)j\f2 (n = 0,±1, ...) 
(see Fig. Etd)). See Fig. Eía)-|SIh) for plots of the spatial and momentum wave-functions of the various approximate 
codewords |0), |1), |+), |— ) in the case (3 = and for the particular vàlues a = 1.8 and r — 1.5. 



B. Intrinsic error probability 

As discussed in Sec.|nJ when approximated codewords are used, one has additional errors (intrinsic errors). In fact, 
due to the presence of the tails of the peaks, the recovery process may lead sometimes to a wrong codeword. The 
recovery in the spatial variable is performed by measuring the operator (i + /3v / 2)(modv / 2)· An intrinsic error in the 
recovery process occurs when, given the state tpo (x) = (x |0), the measurement gives a result within one of the error 
regions: R n = y2 [2n — 3/2 — /3, 2n — 1/2 — fi], n = 0, 1.... The corresponding intrinsic error probability P x $ is equal 

to the one, P x ,i, which we would obtain starting from the state ip± (x) = (ac 1 1) and considering the complementary 
error region. So we simply have 



P X = J2 dx\tp (x)\ 2 . (25) 

„ = 0"'v / 2( 2 «-f-< 3 ) 



It is possible to prové that, under the condition r > 3/2, allowing us to use the asymptotic expansion of the error 
function f °° dte~ l = (2x)~ 1 e~ x [1 — 0(x~ 2 )], we can write the following upper bound for P x : 

P x < lL e -ir+e^/2) [e -o? + Y {vm i + Vm)i] (26) 

jrv 2 . 

v m— 1 

The recovery in the momentum space is instead performed by measuring the operator p(mod7r/\/2). An intrinsic 
error in the recovery process occurs when, given the state ^ + (p), the measurement gives a result within one of the 
error regions fi+ = (tt/ y/2) [2n + 1/2, 2n + 3/2] for integer n, or, given the state ip- (p), the measurement gives a 
result within one of the error regions R~ = (ir/V2) [2n - 1/2, 2n + 1/2] for inte ger n. The corresponding intrinsic 
error probabilities are given by 

+oo „ 

P p-±= E ± d P\^±(p)\ 2 - (27) 



Inserting Eq. (|24|l into Eq. l|27fl we get, after some àlgebra, 



2 



m—0 m>m'=0 



(28) 



where 



V dp exp(-e~ 2r p 2 )(l ± cos V2p) cos[2\/2p(m ~ m')} . (29) 

.6 'Jitè 
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To estimate the quality of the overall encoding procedure provided by this scheme, we have to consider a mean intrinsic 
error probability P e , which is obtained, in general, by averaging over all the possible encoded qubit states. Using the 
above definitions, we have that the mean intrinsic error probability P e satisfies the inequality 

P e < max{P x , P p ,+ ,P p ,-} = P max , (30) 

which defines the maximum intrinsic error probability P max , providing therefore a good characterization of the pro- 
posed encoding scheme. However, in the considered physical configuration, in a large and significant region of param- 
eters, it is P p> + — P p .- = P p 3> P x so that we can take as upper bound for P e simply the quantity P p . 

Let us therefore study the behavior of this upper bound for the intrinsic error probability in the case of system 
parameters corresponding to a typical experimental cavity QED situation. As we have seen, the amplitude a of the 
cavity mode coherent state and the squeezing parameter r the ion CM position are the relevant parameters to study, 
and they should be as large as possible. However, as in the preceding scheme, we have derived our approximatc 
codeword states by making some assumptions and therefore a and r cannot be freely chosen. A first limitation comes 
from the Lamb-Dicke limit. In practice this implies that during all the process of generation of the comb-like state 
the ion has to remain localized within a region smaller than A c so that the linearization of the optical potential is 
vàlid. We can roughly impose \x ± Ax\ < A c /8, where x is the ion mean position and Ax is its mean position spread. 
Using Eqs. I|18|) - (|2Ü|I and l|23|) . we can see that the ion comb-like state <p(x) is substantially confined in the interval 
\/2[2m_ — /3, 2m + — /?], where m± = max(a 2 /2 ± a/y/2,0). Using this fact, we can see that the Lamb-Dickc limit 
implies 

(2m + -/3)|</3g<^, (31) 

which is equivalent to require 



7T VI +4/3-1 
P < ^ = /?ma X , a < - -j= = cw(/3) . (32) 



In order to reach the maximum value of a we choose f3 — /3 max so that 



Vl+7T/(2$)-l _ 

a < a m ax(Pmax) = -J= = «max • (33) 

Therefore we have found an "optimal" value (3 max for /?, and a corresponding upper bound a max for a imposed by the 
adoption of the Lamb-Dicke approximation. Moreover the Lamb-Dicke limit has to be satisfied also by the ion CM 
state at the beginning, and this implies a condition on the squeezing parameter r which is roughly given by r > 3/2. 

Finally, the large detuning approximation imposes another limitation on a. In fact, this condition reads a < S/2go 
and using the condition for having a comb-like state k — 1/2, implying <5 = 2£<7 2 /w a , we get 



a " a — = ^2MüJ! 9 ° kc - (34) 

To summarize, the adopted approximations imply the conditions r > 3/2, a < a'^^ = min(a max , a4 iax ) and — f3 max . 

However, despite these constraints, one can find suitable vàlues, for a and r, such that these bounds are satisfied 
and a low value for the estimated intrinsic error probability is achievable. 

We have done a numerical study of P x , P p .+, Pp.- in the interval < a < 5.5 and for r = 1.5, 2, 3. As anticipated 
above it is P x <C P p — P p .+,P Pi - so that the intrinsic error probability is practically determined by P p . This upper 
bound P p weakly depends on r and, in fact, it has the same behavior for r = 1.5,2,3 , while it has a significant 
dependence upon a. As expected, one has better results for increasing a (see Fig. 0}, even though one cannot take 
too large vàlues for it because of the Lamb-Dicke and large detuning approximations employed here. 

For example, considering the case of a trapped Ca + ion (with Ao — 866nm), if we choose experimentally realistic 
vàlues such as uo a ~ 400KHz and g ~ 3.8MHz, from Eqs. and one gets a'^^ ~ 1 corresponding to 
P — Anax — 1.2. In correspondence of these vàlues for a and /?, and choosing r = 1.5 one has a mean intrinsic 
error probability P e < 9%, as we can see from Fig. 0] These particular vàlues of a and r imply a success probability 
V{X — 0) ~ 26% according to Eq. J2J. Note that, thanks to its exponential-decaying behavior, the value of 
the intrisic error probability can rapidly be improved by accessing stronger ponderomotive interactions (i.e., higher 
go— couplings). 
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FIG. 4: Upper bound P p for the mean error probability P e as a function of a in the interval < a < 5.5 and for r = 1.5 (the 
curve is indistinguishable from those corresponding to r — 2 and r = 3). Numerical data (circles) are fitted by an exponential 
curve (line). 



IV. CONCLUSION 

For any quantum information processing to become a reality the task of providing adequate error correction needs to 
be fulfilled. As the quantum mechanical oscillator is a simply and prevalent model in the study of quantum mechanics, 
it appears to be a natural test bed for such purposes. Thus we have discussed how to embed a qubit in a continuous 
quantum system, so that its redundancy can be used to correct errors which arise from unwanted interactions with the 
environment. In particular we have shown that ponderomotive interaction sumces to this end and we have proposed 
a physical scheme in which a trapped ion within a cavity is considered. 

In this configuration, we arc ablc to show that sufficicntly low vàlues of the intrinsic error probability are reachable 
so that idealized states Q can be effectively engineered with current technology. For the sake of simplicity, throughout 
the paper we have considered protocols conditioned to a given result of a homodyne measurement of the intracavity 
field (the case corresponding to the outcome X = to be specific). However in practical situations one has to consider 
a finite interval of acceptable measurement results. Then, it is possible to see that the success probability is improved 
by enlarging the interval at expenses of an increased error probability. 
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